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Abstract — We consider multi-hop networks comprising Binary 
Symmetric Channels (BSCs). The network carries unicast flows 
for multiple users. The utility of the network is the sum of the 
utilities of the flows, where the utility of each flow is a concave 
function of its throughput. Given that the network capacity is 
shared by the flows, there is a contention for network resources 
like coding rate (at the physical layer), scheduling time (at the 
MAC layer), etc., among the flows. We propose a proportional 
fair transmission scheme that maximises the sum utility of flow 
throughputs subject to the rate and the scheduling constraints. 
This is achieved by jointly optimising the packet coding rates of 
all the flows through the network. 

Index Terms — Binary symmetric channels, code rate selection, 
cross-layer optimisation, network utility maximisation, schedul- 
ing 



I. Introduction 

In a communication network, the network capacity is shared 
by a set of flows. There is a contention for resources among 
the flows, which leads to many interesting problems. One such 
problem, is how to allocate the resources optimally across 
the (competing) flows, when the physical layer is erroneous. 
Specifically, schedule/transmit time for a flow is a resource 
that has to be optimally allocated among the competing 
flows. In this work, we pose a network utility maximisation 
problem subject to scheduling constraints that solve a resource 
allocation problem. 

We consider packet communication over multi-hop net- 
works comprising of Binary Symmetric Channels (BSCs, 
(U). The network consists of a set of C > 1 cells C = 
{1, 2, • • • , C} which define the "interference domains" in the 
network. We allow intra-cell interference {i.e transmissions by 
nodes within the same cell interfere) but assume that there is 
no inter-cell interference. This captures, for example, common 
network architectures where nodes within a given cell use the 
same radio channel while neighbouring cells using orthogonal 
radio channels. Within each cell, any two nodes are within the 
decoding range of each other, and hence, can communicate 
with each other. The cells are interconnected using multi- 
radio bridging nodes to create a multi-hop wireless network. 
A multi-radio bridging node i connecting the set of cells 
B{i) = {ci, .., Cn} C C can be thought of as a set of n single 
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Fig. 1. An illustration of a wireless mesh network with 4 cells. Cells a, 
6, c, and d use orthogonal channels CHi, CH2, CH3, and CH4 respectively. 
Nodes 3, 5, and 6 are bridge nodes. The bridge node 3 (resp. 5 and 6) is 
provided a time slice of each of the channels CHi & CH2 (resp. CH2 & 
CH4 for node 5 and CH2& CH3& CH4 for node 6). Three flows /i,/2, 
and /s are considered. In this example, Cf^ = {a, 6}, Cf^ = {d, b, a}, and 
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radio nodes, one in each cell, interconnected by a high-speed, 
loss-free wired backplane (see Figure [T]). 

Data is transmitted across this multi-hop network as a set J" 
= {1,2,--- ,F}, F > 1 of unicast flows. The route of each 
flow / G J^ is given by Cf = {ci(/), C2 (/),•• • ,Q/(/)}. 
where the source node s{f) G ci(/) and the destination node 
d{f) e C£f{f). We assume loop-free flows (i.e., no two cells 
in Cf are same). Figure [T] illustrates this network setup. A 
scheduler assigns a time slice of duration Tf^c > time 
units to each flow / that flows through cell c, subject to the 
constraint that X1/cgc ^/,c ^ ^c where Tc is the period of the 
schedule in cell c. We consider a periodic scheduling strategy 
in which, in each cell c, service is given to the flows in a 
round robin fashion, and that each flow / in cell c gets a time 
slice of Tf^c units in every schedule. 

The scheduled transmit times for flow / in source cell ci(/) 
define time slots for flow /. We assume that a new information 
packet arrives in each time slot, which allows us to simplify 



the analysis by ignoring queueing. Information packets of each 
flow / at the source node S{f) consist of a block of kf 
symbols. Each packet of flow / is encoded into codewords 
of length Uf = kf/rf symbols, with coding rate < r/ < 1. 
The code employed for encoding is discussed in Section [III 
We require sufficient transmit times at each cell along route Cf 
to allow Uf coded symbols to be transmitted in every schedule 
period. Hence there is no queueing at the cells along the route 
of a flow. 

Channel Model: The channel in cell c for flow / is 
considered to be a binary symmetric channel (BSC) with the 
cross-over probability (i.e., the probability of a bit error) being 
af^c ^ [0, 1]. The corresponding transition probability matrix 
is thus given by 



H/,c(Q^/,c) 



1 - af^c OLf.c 

Thus, the end-to-end channel for flow / is a cascaded 
channel (of £f BSCs), which is a BSC, with the transition 
probability matrix H/(a/) = Hcgc H/,c(<^/,c), the cross- 
over probability of which is given by 
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Since, each transmitted symbol in a packet of a flow can, 
in general, take values from a 2^ = M-ary alphabet, there 
are m channel uses of the BSC for every transmitted symbol. 
Thus, the symbol error probability (for any m > 1) is given 
by /3f = 1 — (1 — Of/ )^. Let the Bernoulli random variable 
Ef [i] indicate the end-to-end error of the ith coded symbol 
at the destination in a code word of flow /. Note that Ef[i]s 
are independent and identically distributed (i.i.d.), and that 
P{Ef\i] = 1} = f3f = l-P{Ef\i] = 0}. In the channel model 
described, the channel processes across time are independent 
copies of the BSCs. This is realised in a wireless network 
by means of an interleaver of sufficient depth (after the 
channel encoder), which interleaves the encoded symbols. The 
interleaved symbols see a fading channel (which is modelled as 
a channel with memory, e.g., a Gilbert-Elliot channel ||2l), but 
the de-interleaver (before the channel decoder) brings back the 
original sequence of the encoded symbols, but interleaves the 
channel fades, the combined effect of which can be modelled 
as independent channel processes across time. In another 
work 1 3 1, we model the fading channel as a packet erasure 
channel (or a block fading channel), and obtain the optimal 
transmission strategy, which includes optimal interleaving of 
bits across schedules and the optimal coding rates. 

Letting e/ (r/) denote the error probability that a packet fails 
to be decoded, the expected number of information symbols 
successfully received is Sf{rf) = /c/(l — e/(r/)). Other things 
being equal, one expects that decreasing r/ (i.e., increasing 
the number of redundant symbols Uf — kf) decreases error 
probability e/, and so increases Sf. However, since the 
network capacity is limited, and is shared by multiple flows, 
increasing the coded packet size n/^ of flow /i generally 
requires decreasing the packet size Uf^ for some other flow 
/2. That is, increasing Sf^ comes at the cost of decreasing 



Sf^ . We are interested in understanding this trade-off, and in 
analysing the optimal fair allocation of coding rates amongst 
users/flows. 

Contributions: Our main contribution is the analysis of 
fairness in the allocation of coding rates between users/flows 
competing for limited network capacity. In particular, we pose 
a resource allocation problem in the utility-fair framework, 
and propose a scheme for obtaining the proportional fair 
allocation of coding rates, i.e. the allocation of coding rates 
that maximises "^ f ^jr log Sf{rf) subject to network capacity 
constraints (or scheduling constraints). Specifically, at the 
physical layer, the (channel) coding rate of a flow can be 
lowered (to alleviate its channel errors) only at the expense of 
increasing the coding rates of other flows. Also, at the network 
layer, the length of schedules of each flow should be chosen in 
such a way that it maximises the network utility. Interestingly, 
we show in our problem formulation that the coding rate and 
the scheduling are tightly coupled. Also, we show that for a log 
(network) utility function (which typically gives proportional 
fair allocation of resources) the optimum rate allocation (in 
general) gives unequal air-times which is quite different from 
the previously known result of proportional fair allocation 
being the same as that of equal air-time allocation (H). This 
problem, which we show in Section [IIIl requires solving a 
non-convex optimisation problem. Our work differs from the 
previous work on network utility maximisation (see |5] and 
the references therein) in the following manner. To the best of 
our knowledge, this is the first work that computes the optimal 
coding rate for a given scheduling (or capacity) constraints in 
the utility-optimal framework. 

The rest of the paper is organised as follows. In Section [III 
we obtain a measure for the end-to-end packet decoding error, 
and describe the throughput of the network. In Section |III1 we 
formulate a network utility maximisation problem subject to 
constraints on the transmission schedule lengths. We obtain 
the optimum coding rates for each flow in the network in 
Section |lVl In Section |Vl we provide some simple examples 
to illustrate our results. The proofs of various Lemmas are 
omitted due to lack of space. 

II. Packet Error Probability 

We recall that each transmitted symbol of flow / reaches 
the destination node erroneously with probability /3/. Hence, 
to recover the information packets, we employ a block code 
at the source nodes (a convolutional code with zero-padding 
is also a block code). Since an (n^k^d) code can correct up 
to [^^J errors, we are interested in employing a code with a 
large distance d. Thus, a natural choice is the class of (linear) 
maximum-distance separable (MDS) codes. MDS codes of 
rate k/n have the property that it achieves the Singleton bound 
(L6J), 

d^n-k^l, (1) 

i.e., the minimum distance between any two codewords d, in 
an MDS code is n — k -\- 1. Thus, the maximum number of 
errors that an MDS code can correct is [^^J = [^^-^J • It is 
well known that in the case of binary signalling, only trivial 



MDS codes exist. Hence, in this paper, we consider M = 2^- 
ary alphabet, where m > 1. Examples for MDS codes in the 
case of non-binary alphabets include Reed-Solomon codes 
(IH), and MDS-convolutional codes (|7]). In f?], the authors 
show the existence of MDS-convolutional codes for any code 
rate. We note here that Reed-Solomon codes can also correct 
burst errors, and hence, is more suitable for wireless networks 
(which does not employ an interleaver). 



A. Network Constraints on Coding Rate 

Based on the modulation and the bandwidth available at 
each cell c, a flow /, which passes through it, can obtain a 
maximum feasible physical (PHY) rate of transmission in bits 
per second that the cell c can support. Let Wf^c be the PHY 
rate of transmission of flow / in cell c. For each transmitted 
packet of flow /, each cell c e Cf along its route must allocate 
at least ^^ units of time to transmit the packet (or encoded 
block) where we recall that n/ is the length of the code word. 
Let Tc := {/ G J^ : c G C/} be the set of flows that are routed 
through cell c. We recall that the transmissions in any cell c 
are scheduled in a TDMA fashion, and hence, the total time 
required for transmitting packets for all flows in cell c is given 
by y'fcTr ^^- Since, for cell c, the transmission schedule 
interval is Tc units of time, the coding rates r/ must satisfy 



the schedulability constraint ^ 



feJ^c rfZ 



< T 



B. Error Probability - Upper bound 

The symbol errors Ef[l]^Ef[2]^- - ■ ^Ef[nf] are i.i.d. 
Bernoulli random variables, and hence, the probability of 
a codeword (or encoded packet) being decoded incorrectly 
is given by P < Yl^^i Ef[i] > ^^~ ^ \. We observe that 

Yll^=i^f[^] is ^ binomial random variable, and hence, the 
probability of decoding error can be computed exactly. How- 
ever, the exact probability of error is not tractable for further 
optimisation as the probability of error, which is a function of 
the coding rate, is neither concave nor convex. Hence, we pose 
the problem based on the upper bound on the error probability 
So, we obtain an upper bound and a lower bound for the error 
probability. We show that the bounds are tight, and hence, the 
problem of network utility maximisation can be posed based 
on the lower bound on the error probability. 

Lemma 1. An upper bound for the end-to-end probability of 
a packet decoding error for flow f is bounded by the following. 



e/= P<^^^/W> 



nf-kf 



i=l 






(2) 



where Of > is the Chernoff-bound parameter and the 
function Iz{x;6) := 6x — ln(E [e^^]) is called the rate 
function in large deviations theory. 



C Error Probability - Lower bound 

Lemma 2. The end-to-end probability of a packet decoding 
error for flow f is at least as large as 

1^1 ...^( kf 



e/ 



> 



1-/3/ 



• exp 



exp 
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l-2xf 



D{B{xf)\\B{Pf)) 



(3) 



where B{x) is the Bernoulli distribution with parameter x, 
H(V) is the entropy of probability mass function (pmf) V, 
and D{V\\Q) is the information divergence between the pmfs 
V and Q. 

From the lower and the upper bounds for the probability of 
packet decoding error, and for the optimal 6>j (see Eqn. ([T5]) 
in Section HVl) , we see that the exponent of the lower bound 
is the same as that of the upper bound (Eqn. (O) with a pre- 
factor. This motivates us to work with the lower bound e/ as 
a candidate to compute the utility of flow /, which is given 
byln(fc/(l-e/)). 

We recall that Ef[l] is a Bernoulli random variable which 
takes 1 with probability /3/, and with probability 1 — /3/. 

Thus Ie,w (^;^/) = Of (^) - In (1 - /?/ + /?/e^^l. 
Let Xf : 
Eqn. ©, 



-^^. Note that ^ x/ < |. Therefore, from 



CfiOf^Xf) := exp (- J^^ i^f^f - In (l - /3/ + pfc'f)] 

(4) 



in. Network Utility Maximisation 

We are interested in maximising the utility of the network 
which is defined as the sum utility of flow throughputs. 
We consider the log of throughput as the candidate for the 
utility function being motivated by the desirable properties 
like proportional fairness that it possesses. 

We define the following notations: Chernoff-bound parame- 
ters 6 := [OfjfeT, code rates r := [r/J/^j^, and x parameters 
X := [xf]feT (where we recall that Xf = {1 — r/)/2). We 
define the network utility as 

U{0,x) := ^ln(fej(l-e/(^/,x/))) 

= ^ In (kf) + ^ In (1 - ef{Of, Xf)) . (5) 



/e-F 



/e-F 



The problem is to obtain the optimum coding rate parameter 
X* and the optimum Chernoff-bound parameter ^*, which 
maximises the network utility. Since, kf , the size of informa- 
tion packets of each flow / is given, maximising the network 
utility is equivalent to maximising 

U{0,x) := ^\n{l -Cf {Of, Xf)). (6) 

Thus, we define the following problem 



PI: 



max [/(0,x)=^ln(l-e/(^/,x/)) 
kf 



subject to XI Ti 



f:ceCf 

Xf < A/ 
Xf > A/ 



(1 - 2xf)wf,, 



<Tc. 



\/ceC (7) 



(8) 
We note that the Eqn. ^ enforces the network capacity 

(or the network schedulabiHty) constraint. The objective 

function U{6^x) is separable in {Of^Xf) pair for each flow 

/. Importantly, the component of utility function for each 

flow / given by In (1 — e/(^/, x/)) is not jointly concave in 

(Of^Xf). However, In (1 — e/(6>/,x/)) is concave in Of (for 

any Xf), and in Xf (for any Of). Hence, the network utility 

maximisation problem PI is not in the standard convex 

optimisation framework. Instead, we pose the following 

problem. 



P2: 



max max N In (1 — e/(^/, x/)) 



(9) 



/e-F 



subject to 



kf 



(1 - 2xf 

Of>0, 
Xf <Xf 
Xf >^f 



^^f,c 



< Tc, \/ceC 

V/G^ 



(10) 
In general, the solution to P2 need not be the same as the 
solution to PI. However, in our problem, we show that P2 
achieves the solution of PI. 

Lemma 3. . For a function f : y x Z ^ R that is concave 
in y and in z, but not jointly in {y^ z), the solution to the joint 
optimisation problem for convex sets y and Z 



is the same as 



max f (v, z) 



maxmax/(?/, z), 

zez yey 



(11) 
(12) 



if f{y*{z)^ z) is a concave function of z, where for each z ^ Z, 
y*{z) :=argmax/(7/,z). 

We note that for each x/, the probability of error ef{Of^Xf) 
is convex in Of, and hence, ln(l — e/) is concave in Of. Thus, 
we first solve for the optimum Chernoff bound parameter 0^ 
which we describe in Section IIV-A[ After having solved for 
the optimum ^*, we show in Section HV-Bl that U{0*(x), x) is 
a concave function of x. Hence, from Lemma [Sj the solution 
to problem (P2) (the maximisation problem that separately 
obtains the optimum ^* and optimum a?*) is globally optimum. 
We study the rate optimisation problem that obtains cc* in 
Section [I^Cl 



IV. Utility Optimum Rate Allocation 
A. Optimal 6>* 

Consider the following optimisation problem, for any given 

xe[\f,\fY. 

max ^ In (1 - ef{Of,Xf)) (13) 

subject to Of > 0, yfeJ^ 

We note that the objective function is separable in OfS, and that 
Cf is convex in Of. Hence, the problem defined in Eqn. ([13]), 
is a concave maximisation problem. We recall that 

kf 



ef{Of,Xf)=exp 



l-2xf 



[OfXf -In {1-Pf + Pfe'f)] 
(14) 



The partial derivative of ej with respect to 0f is given by 



def_ 

Observe that 
if, for Of 



e/ 



y 



l-2xf 



Xf 



Pfe' 



1 - /3/ + Pfe^f 



111 



0,Xf - J 



9j is an increasing function of Of. Thus, 

7/ = u, Xf - i_s\b < or Xf < Pf (equivalently, 
r/ > 1 — 2/3/), the derivative is positive for all Of > 0, or 
6/ is an increasing function of Of. Hence, for Xf < Pf, the 
optimum 0% is arbitrarily close to which yields Cf arbitrarily 
close to 1. Thus, for error recovery, for any end-to-end error 
probability /3/, the coding rate should be smaller than 1 — 2/3/, 
in which case, we obtain the optimal 0*n by equating the partial 
derivative of e/ with respect to Of to zero. 



I.e., 



or, 
or, 

The probability of error for a given Xf and 0*f{xf) is then 
given by 



/J/e'^ _ ^_ 






-i-(S; 



ef{0},Xf) 



exp 



exp 



kf 



-2xf [ 
kf 



Xf\n{^ 



(l-a;/)ln 



^f 



1-/3/ 



l-2xf 



Dmxf)\\B{Pf) 



(15) 



B. A convex optimisation framework to obtain optimal x*^ 

If ln(l — ef{0*Axf), Xf)) is a concave function of Xf, then 
one can obtain the optimum x*f using convex optimisation 
framework. To show the concavity of ln(l — ef{6Uxf), Xf)), 
it is sufficient to show that ef{6j{xf),Xf) is convex in Xf. 

Define A/ := In (f4^^] . Note that 



def 
dxf 

dxi 



-ef 



Pt{i-l3t)J- 
kfAf 



e/ 



{l-2xf)^ 

kf ' 

(l-2x/)2J 



kf 



(1 - 2xf) 



A2 
2^V 



4A/ 



2xf 



l-2xf Xf{l-Xf) 



ef{OUxf)^Xf) is convex if 



kf 



(1 - 2xfy 



-Ai> 



4Ai 



2xi 



l-2xf Xf{l-Xf)' 



or, 



4(1 -2x/) , (1 - 2xf)' 



< 



kf 



Af Xf{l-Xf)A} 

Since, we consider Xf ^ Aj, where Aj = /3f 



^f 



for 



some arbitrarily small e/ > 0, we have ^ ^ Kq where 

I/Kq := In f ~^^ nj^^ ), and hence, a sufficient condition for 
the convexity of e/ (and hence, the concavity of ln(l — e/)) 
is 



4(1 - 2xf) 



K, 



2(l-2a;/f 



< 



kf 



(16) 



A/ ''°xf{l-xf) 

The above condition is a convex function of x/, and we 
include this as a constraint in the problem formulation. Thus, 
ef{0*^{xf)^Xf) is convex in x/, and hence, we obtain the 
optimal Xj using convex optimisation method. Also, from 
Lemma [3l the optimal coding rate r j = 1 — 2x% is unique 
and globally optimum. 

The minimum kf required to ensure convexity of 
ef{OUxf)^Xf) is computed numerically, and is tabulated 
below. 

TABLE I 
Minimum kf that ensures convexity of ef{6%{xf), Xf) 



^f 


minimum kf required 


0.1 
0.01 
0.001 
0.0001 


6 
10 
33 
164 



From the above table, we see that the minimum packet size 
required to ensure convexity is very small, and in practice, 
the packet size kf is much larger than the minimum size re- 
quired. Hence, for all practical purposes, the optimal code rate 
problem is a convex problem. More importantly, the constraint 
given by Eqn. ([T6b is not an active constraint. However, for 
the sake of completeness, we include this constraint in the 
problem definition below. 

C Optimal Coding Rate r 

In this subsection, we obtain the optimal coding rate using 
the optimal Chernoff-bound parameter vector ^*, by solving 
the following network utility maximisation problem 



max^ln(l-e/((9},x/)) 
Xf < A/ 

Xf >\j 



(17) 

^ceC 

4(l-2.,) ,^, (l-2.,)3 ^^^ V/.^ 



A 



/ 



'^/(l-^/) 



(18) 



The objective function is separable and concave, and hence, 
can be solved using Lagrangian relaxation method. Also, the 
constraint represented by Eqn. (fTSl) is not an active constraint, 
and hence, there is no Lagrangian cost to this constraint. We 
note here that the coding rate should be such that kf/{l — 2xf) 
is an integer, and hence, obtaining x j is a discrete optimisation 
problem. This is, in general, an NP hard problem. Hence, we 
relax this constraint, and allow Xf to take any real value in 
[A J, A/]. The Lagrangian function for the optimal rate problem 
is thus 



L{x,p,u,v) 



i:,„,-,.>,.,)-g.(E^ 



2xf)wf^c 



Tr, 



Y^ Uf {xf - A^) - ^ Vf {xf - Xf) 

feT feT 



Applying KKT condition, ^^ U*= 0, we have 

Pc '^kf 

ceCf 

2kf I ^^ p 



-1 9ef 

1 - ef dxf '"^ ;^^ Wf^c (1 - 2xp2 



EPc ^kf 



*\2 



(1 - 2x}) 



IZ ^7^ + ^/ - ^/ 



.ceCf 



Wf, 



e/ kfA) 



1-ef (1 - 2xp2 (1 _ 2x}r \^^ Wf, 



i^^J-MS 



Pc 1 , JVf - ^/)(1 - ^X})' 



e/ 



<ceC, 



Wf,c 



y 



*\2 



(vf - Uf)il - 2x}) 



kf 

where Xf := 2 (^eec, ,^) and A} := In (g{^). If 
the optimal x*. is either A^ or A/, then it is unique. If x*. e 
(A/, A/), then Uf = Vf = 0, and in this case (which is the 
most interesting case, and we consider only this case for the 
rest of the paper), we have 



e/ 



1-ef 



.A} = Xf 



A 



e/ 



/ 



A/+A} 
A/ 



(19) 



^^M-r^^(^(^/)ii^(^^»)^A,+A^ 



^^L_D{B{x*f)\\B{l3f)) = \n 



A/ + A} 
A/ 



(20) 



In the above equation, both the LHS and the RHS are 
increasing in x*f. Also, LHS is a strictly convex (increasing) 
function and RHS is a strictly concave (increasing) function 
of x*f. Hence, they intersect at exactly one point in the region 
(/3/, 0.5] which is the optimal x*. for a given Lagrangian price 
vector p. 



D. Sub-gradient Approach to Compute optimal p* 

In this section, we discuss the procedure to obtain the 
optimal shadow costs or the Lagrange variables p* . The dual 
problem for the primal problem defined in Eqn. (fTTl) is given 
by 

min D(p), 

where the dual function D{p) is given by 
D{p) 



Y: ln(l - ef{x){p))) + Y.PC I T. - Y: TYZr^ 



(21) 



(p))^/,< 



(22) 



In the above equation, ef{xf) denotes ef{OUxf)^Xf). Since 
the dual function (of a primal problem) is convex, D is convex 
in p. Hence, we use a sub-gradient method to obtain the 
optimum p* . From Eqn. ([2T1) . for any x, 

Dip) > Y ln(l - cfi.f)) + Yp^ U - E n-2t'w. ) ' 

and in particular, the dual function D{p) is greater than that 

for X = x*Ap), i.e., 

D{p) 
Thus, a sub-gradient of !)(•) at any p is given by the vector 

k, 



^'-^n 



feTc 



(1 - 2x}ip))wf,c 



(24) 



cec 



We obtain an iterative algorithm based on sub-gradient method 
that yields p*, with p{i) being the Lagrangians at the ith 
iteration. 



Pcii^l) 



Pc{i) - 7 • \Tc 



k 



f 



2^{i-2x}{pmwf,c 



feT. 



where 7 > is a sufficiently small stepsize, and [f{x)\^ := 
max{/(x), 0} ensures that the Lagrange multiplier never goes 
negative. Note that the Lagrangian updates can be locally done, 
as each cell c is required to know only the rates x*Ap{i)) of 
flows j ^ Tc- Thus, at the beginning of each iteration i, the 
flows choose their coding rates to 1 — 2x^(p(z)), and each cell 
computes its cost based on the rates of flows through it. The 
updated costs along the route of each flow are then fed back 
to the source node to compute the rate for the next iteration. 



The Lagrange multiplier pc can be viewed as the cost of 
transmitting traffic through cell c. The amount of service time 
that is available is given by A = Tc — V .-, ^ >n/ ^^/•^^^ . 

to -^ c ^ (l-2x*(p(2)))'u;/,c 

When A is positive and large, then the Lagrangian cost Pc 
decreases rapidly (because D is convex), and when A is 
negative, then the Lagrangian cost Pc increases rapidly to 
make A > 0. We note that the increase or decrease of pc 
between successive iterations is proportional to A, the amount 
of service time available. Thus, the sub-gradient procedure 
provides a dynamic control scheme to balance the network 
load. 

We explore the properties of the optimum rate parameter 
X J in Section IIV-EI In Section |Vl we provide some examples 
that illustrate the optimum utility-fair resource allocation. 

E. Properties ofx""^ 

We are interested in studying the behaviour of the optimum 
coding rate r^ = 1 — 2x%, when the PHY rate Wf^c and 
the packet size kf increases such that kf/wf^c is always a 
constant. 

Lemma 4. r^ = 1 — 2x*Jkf) is an increasing function of kf 
(with the PHY rate Wf^c being proportional to kf). 

Lemma m is quite intuitive. For any given channel error /3f, 
as the block (or packet) length increases, it is optimum to go 
for a high rate code. In other words, it is optimum for a flow 
to use as much scheduling time as possible (i.e., use a large 
block length kf, and hence, use a high rate code); however, 
the resources are shared among multiple flows, and hence, we 
ask the following question: ''what is the optimum share of the 
scheduling time'' that each flow should have. Interestingly, in 
our problem formulation, the optimum code rate parameter Xj 
also solves this optimum scheduling times for each flows. 

It is interesting to ask the question of how large the packet 
sizes kf be for optimum resource allocation, and Lemma |4] 
provides a hint to the solution. From Lemma (H we understand 
the following: if there are two flows /i,/2, through a cell c 
(seeing the same channel conditions, i.e., /3f^ = f3f^) with 
^/i,c > ^/2,c then it is optimum for flow /i to use a large 
packet size kf^ and flow /2 to use a small packet size kf^. 
The optimum schedule length will be to allocate less schedule 
time to flow /i and more schedule time to flow /2. 

In the asymptotic case when Wf^c <^nd kf grows to 00 (and kf 
grows linearly with Wf^c, we see from Eqn. (|2Q|) that the error 
exponent also goes to 00 (as 1 — 2xf > 0), and hence, e/ -^ 
0. In this case, we see that the optimum rate can approach 
arbitrarily close to 1 — 2/3j. Thus, for any kf and t^/,c, the 
optimum coding rate r^ < 1 — 2/3^ 

Previous studies on optimum resource allocation establish 
that the proportional fair allocation is the same as equal 
air-time allocation ((H). But, in this problem, we see an 
interesting phenomenon that is unusual of a proportional-fair 
resource allocation. 

Lemma 5. The optimum rate allocation x* (or equivalently 
r*j is not equivalent to equal air-time allocation which is 





Fig. 2. Cells with equal traffic load 



Fig. 3. Cells with unequal traffic load 



typically the solution of a proportional-fair (or In utility) 
allocation. 

In particular, we see that the flows that see a better channel 
get less air-times than the flows that see a worse channel. 
This phenomenon is evident in the case of infinitely long code 
words; with other parameters being same, the air-times of 

, and hence, flows 



■2/3/,c 



flows in a cell c are proportional to — 
with small (3 get less air-times. 

V. Examples 

In this Section, we analyse some simple networks based on 
the utility optimum solution that we obtained. In particular, 
we analyse the so-called parking-lot topology often used to 
explore fairness issues. It is to be noted that the parking-lot 
topology is a simple case of a line network, and the results of 
this section extends in a simple way to a linear network. 

A. Example 1: Two cells with equal traffic load 

We begin by considering the example shown in Figure [21 
consisting of two cells a and h having three nodes 1, 2, and 
3. Each cell has the same symbol error probability P and the 
schedule length T. There are three flows /i, /2, and /a, with 
two of the flows /i and /s having one-hop routes C/^ = {^} 
and C/3 = {a}, and one flow /2 having a two-hop route 
C/2 = {a, h). Each flow has the same information packet size 
k and PHY transmit rate, i.e. Wf^c = ^• 

The end-to-end packet error probability experienced by the 
two-hop flow /2 is greater than that experienced by the one 
hop flows /i and /a, since each hop has the same fixed error 
probability. Hence, we need to assign a lower coding rate rf^ 
to flow /2 than to flows /i and fs in order to obtain the 
same error probability (after decoding) across flows. However, 
when operating at the boundary of the network capacity region 
(thereby maximising throughput), decreasing the coding rate 
r/2 of the two-hop flow /2 requires that the coding rate of 
both one-hop flows /i and fs be increased in order to remain 
within the available network capacity. In this sense, allocating 
coding rate to the two-hop flow /2 imposes a greater marginal 
cost on the network (in terms of the sum-utility) than the one- 
hop flows, and we expect that a fair allocation will therefore 
assign higher coding rate to the two-hop flow /2 . The solution 
optimising this trade-off in a proportional fair manner can be 
understood using the analysis in the previous section. 

In this example, both the cells are equally loaded and, by 
symmetry, the Lagrange multipliers Pa = Pb- Hence, A/^ = 



A/2 

find from Eqn. ([T9l) that 



-J2. = Xr Note that x% < x% and A 



/i 



V2 



< Aj^. Hence, 



we 



e/2 



A/, A/2 



•A* 



A/2 A/, + A*^ 
< 1. 



B. Example 2: Two cells with unequal traffic load 

We consider the same network as in the previous example, 
but now with only the flows /i and /2 (i.e., the flow /s is 
not present, see Figure [3]) in the network. In this example, cell 
h carries two flows while cell a carries only one flow. The 
encoding rate constraints are given by 

1 wT 

— < — — , (from cell a), 

r/2 k 

11 wT 
\ < — — , (from cell h). 

Since, both r/^ and Vf^ are at most 1, it is clear that at 
the optimum point, the rate constraint of cell a is not tight 
while the constraint of cell h is tight. Thus, the shadow prices 
(Lagrange multipliers) Pa = and p6 > 0. That is, at the 
first hop the cell is not operating at capacity, and so the 
"price" for using this cell is zero. In this example, A/^ = A/2, 
and hence, from Eqn. ([T9l ), we deduce that for low channel 
errors, e/^ ^ e/2 . This allocation make sense intuitively since 
although flow /2 crosses two hops, it is only constrained at the 
second hop and so it is natural to share the available capacity 
of this second hop approximately equally between the flows. 



VI. Conclusions 

In this paper, we posed a utility fair problem that yields 
the optimum coding across flows in a capacity constrained 
network. We showed that the problem is highly non-convex. 
However, we provided some simple conditions under which 
the global network utility optimisation problem can be solved. 
We obtained the optimum coding rate, and analysed some of its 
properties. We also analysed some simple networks based on 
the utility optimum framework we proposed. To the best of our 
knowledge, this is the first work on cross-layer optimisation 
that studies optimum coding across flows which are competing 
for network resources. 
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